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Abstract 


\ 

A  Horizontal  Line  Array  (HLA)  is  subject  to  left-right  bearing  ambiguity  and  also  to 
bearing  bias  caused  by  non-zero  vertical  arrival  angles.  However,  if  a  second  observation  is 
made  with  the  array  rotated  to  a  different  orientation,  the  azimuthal  and  vertical  arrival  angles 
may  be  estimated.  For  two  observations  where  the  source-array  geometries  are  static,  a  closed- 
form  solution  is  available.  For  the  solution  error,  a  closed-form  approximation  is  also  available, 
but  only  if  the  observation  errors  are  small.  This  report  deals  with  the  case  of  moderate  to  large 
observation  errors  and  presents  error  bounds  on  the  estimated  azimuthal  and  vertical  arrival 
angles.  As  well,  several  examples  are  included.  l  \  <  ,  /  ■  , 


Sommaire 

Un  reseau  lineaire  horizontal  (HLA)  peut  comporter  un  ambiguite  vers  la  gauche  ou  la 
droite  ainsi  qu’un  biais  de  gisement  en  raison  des  angles  d’arrivee  verticaux  non  nuls.  Toutefois, 
si  on  effectue  une  seconde  observation  apres  avoir  amene  le  reseau  a  une  autre  orientation,  on 
peut  estimer  les  angles  d’arrivee  azimutaux  et  verticaux.  Si  l’on  dispose  de  deux  observations 
dans  lesquelles  la  geometrie  source-reseau  est  statique,  on  peut  obtenir  une  solution  de  forme 
close.  On  peut  aussi  obtenir  une  solution  de  forme  close  pour  la  fonction  d’erreur  de  la  solution, 
mais  seulement  si  les  erreurs  d’observation  sont  d’un  tres  faible  niveau.  Le  rapport  traite  des 
cas  d’erreurs  d’observation  de  niveaux  moyens  a  tres  eleves  et  presente  les  limites  d’erreur  sur 
les  angles  d’arrivee  azimutaux  et  verticaux  estimes.  II  contient  aussi  plusieurs  exempies. 
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1  Introduction 


A  Horizontal  Line  Array  (HLA)  is  a  one  dimensional  array  of  sensors  lying  in  a  horizontal 
plane.  Such  an  array  may  be  used  to  detect  and  estimate  the  arrival  direction  of  a  signal.  The 
use  of  an  arrival  angle  as  bearing  estimate  causes  two  major  problems:  the  arrival  angle  does  not 
contain  any  information  with  regard  to  which  side  of  the  array  the  signal  is  arriving  from(this 
is  the  left-right  ambiguity  problem);  and  if  an  arrival  angle  has  a  non-zero  vertical  component, 
the  bearing  will  be  biased  toward  broadside.  For  arrays  which  can  be  rotated  in  the  horizontal 
plane,  a  technique  has  been  developed  to  resolve  these  problemsjl].  Should  there  be  errors  in 
the  observed  arrival  angles,  the  resulting  estimates  will  show  a  corresponding  error.  The  main 
thrust  of  this  memorandum  is  to  investigate  the  effects  such  errors  have  on  the  estimates  of 
azimuthal  and  vertical  arrival  angles. 

The  following  two  sections  discuss  the  bearing  bias  problem  and  the  prescribed  technique  to 
solve  it.  Section  2  shows  that  the  bias  is  more  pronounced  towards  the  end  fire  directions(signals 
arriving  from  either  end  of  the  array)  compared  to  the  broadside  direction.  The  technique 
described  in  section  3  allows  the  true  azimuthal  and  vertical  arrival  angles  to  be  estimated  if 
the  array  is  rotated  in  the  horizontal  plane  to  a  different  orientation  and  a  second  arrival  angle 
observation  is  made. 

The  main  contribution  of  this  memorandum  is  found  in  section  4,  which  deals  with  error 
bounds  on  the  estimated  azimuthal  and  vertical  arrival  angles.  Given  arbitrary  error  bounds 
on  the  arrival  and  rotation  angles,  analytic  bounds  on  the  estimates  of  azimuthal  and  vertical 
arrival  angles  are  derived. 

In  Section  5,  the  analytic  error  bounds  on  the  estimated  azimuthal  and  vertical  arrival 
angles  are  used  to  validate  the  case  of  small  observation  errors  in  reference  [1]  and  to  illustrate 
the  discrepancies  when  the  observation  errors  violate  the  small  error  restriction. 
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2  The  HLA  Bearing  Bias  Problem 


Figure  2.1  illustrates  the  geometry  of  a  plane  wave  signal  arriving  at  a  Horizontal  Line 
Array(HLA)  lying  in  the  plane  <f>  =  0  where  <f>  is  the  vertical  arrival  angle.  The  azimuthal 
arrival  angle  6  of  the  signal  is  measured  clockwise  from  the  HLA  forward  axis  direction.  For 
such  an  array,  two  problems  arise:  the  first  is  the  so-called  left-right  bearing  ambiguity;  and 
the  second  is  the  bearing  bias.  These  two  problems  are  clarified  below. 


Figure  2.1:  Arrival  Angle  Geometry 


The  HLA  signal  processing  provides  an  estimate  of  /?,  the  angle  of  incidence  between  the 
array  axis  and  the  signal  vector.  As  a  function  of  the  azimuthal  angle  and  the  vertical  angle, 
the  angle  /?  is  given  by 

cos  f3  =  cos  9  cos  4>.  (2.1) 

From  Equation  (2.1),  the  angle  of  incidence  is  given  by 

fi  =  ±cos-1  (cos  0  cos  ^)  (2.2) 

which  has  a  left-right  ambiguity.  Thus,  if  <p  is  restricted  to  lie  between  0°  and  90°,  the  bearing 
9  also  has  the  right-left  ambiguity. 
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If  the  signal  is  in  the  plane  <f>  =  0°,  then  cos  0  =  cos  0.  However,  if  <f>  ^  0°,  cos  0  is  a  biased 
estimator  of  cos0.  The  effect  of  non-zero  <p  is  to  make  |  cos 0\  smaller  than  |  cos  9 j  so  that  9  is 
biased  towards  the  broadside  direction  (where  cos  0  =  0). 

To  illustrate  the  seriousness  of  the  bias  B  which  is  defined  by 

B  =  0-9 

=  cos-1(cos0cos<£)  -  9,  (2.3) 

Figure  2.2  shows  values  of  B  as  a  function  of  0  for  several  values  of  <f>  between  0°  and  50°.  From 
the  curves  in  Figure  2.2,  it  is  evident  that  the  vertical  arrival  angle  can  have  a  very  significant 
effect  and  that  the  maximum  bias  occurs  in  the  endfire  directions  0  =  0°  and  0  =  180°. 


Figure  2.2:  HLA  Bearing  Bias  as  a  function  of  azimuthal  angle  0,  for  several  values  of  vertical 
arrival  angle  <t>. 
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3  Estimation  of  Azimuthal  and  Vertical 
Angles 


3.1  Method  of  Solution 

The  true  values  of  the  azimuthal  and  vertical  arrival  angles  may  be  estimated  if  a  sec¬ 
ond  observation  is  made  with  the  Horizont'l  Line  Array  rotated  in  the  horizontal  plane  to  a 
different  orientation.  For  two  such  observations,  the  estimation  problem  can  be  modelled  by 
two  nonlinear  equations  in  two  unknowns.  A  closed  form  solution  for  this  case  is  possible.  In 
this  section,  the  assumptions  and  restrictions  on  the  model  are  presented,  the  analytic  model 
is  described,  and  the  analytic  solution  of  the  azimuthal  and  vertical  angles  is  given. 


3.2  Assumptions  and  Restrictions 

To  find  a  computable  closed-form  solution,  the  following  simplifying  assumptions  are  made: 
1.  The  array  is  assumed  to  be  perfectly  stable:  that  is,  the  errors  due  to  tilting  and  stretching 
of  the  array  are  ignored.  2.  The  effects  of  motion  between  observations  is  ignored:  that  is, 
the  source  and  receiver  are  assumed  to  be  stationary.  3.  The  propagation  path  is  assumed  to 
be  the  same  for  both  observations:  that  is,  the  vertical  arrival  angle  <j>  is  assumed  to  remain 
constant. 

Regarding  the  above  assumptions,  the  implications  are  as  follows:  1.  Array  tilting  may  be 
included  by  enlarging  the  error  interval  for  the  array  observations.  2.  The  error  intervals  for 
the  input  data  may  be  enlarged  to  compensate  for  the  relative  motion  between  the  source  and 
array.  3.  The  assumption  of  constant  ip  restricts  the  usefulness  of  the  method  to  propagation 
paths  with  stable  vertical  arrival  angles. 


3.3  The  Analytic  Model 

Let  the  Horizontal  Line  Array  be  rotated  in  the  horizontal  plane  through  an  angle  AO 
measured  in  a  clockwise  direction  and  let  the  observed  angles  of  incidence  before  and  after  the 
turn  be  fa  and  fa  respectively.  Let  &i  and  9 2  be  the  corresponding  azimuthal  angles  before 
and  after  the  turn  and  <f>i  and  fa  the  corresponding  vertical  arrival  angles.  Then  it  is  possible 
to  describe  the  model  by  two  equations  and  four  unknowns:  namely,  by 


(3.2) 


COS  @2  =  cos  9 2  COS  <f> 2 .  (3.2) 

The  four  unknowns(0i,  02 ,  4>i,  and  $2)  may  be  reduced  to  two(0  and  <f>)  by  assuming 
constant  vertical  arrival  paths  (<?t  =  <f>2  =  <f>)  and  constant  horizontal  arrival  paths  (#i  =  8, 
and  82  =  6-  A 9)  where  A 6  is  the  known  rotation  angle.  Under  these  assumptions  equations 
(3.1)  and  (3.2)  become 

cos  Pi  =  cos  6  cos  <f>  (3-3) 

and 

cos  P2  =  cos  (0  —  A0)  cos  <f>.  (3-4) 

The  vertical  angle  <p  may  be  restricted  to  the  interval  [0,90]  while  the  azimuthal  angle  8 
may  be  restricted  to  (-180,180].  Both  of  the  incident  angles  pi  and  p 2  may  be  restricted  to 
[0.180],  The  rotation  angle  A9  may  be  restricted  to  (-180,180)  excluding  the  point  A9  =  0. 

3.4  Analytic  Solution 

The  above  system  of  equations  may  be  solved  by  expanding  equation  (3.4)  and  combining 
it  with  equation  (3.3)  [l].  The  solutions  for  <f>  and  8  are  given  by 


4>{Pi,P2,  A0)  =  cos  1\j(^ 


cos  pi  cos  A  8  -  cos  P2 
sin  A  8 


^  +  cos2  Pi 


ana 

•(ft,*. A«)  =  »jn(sm#(A,ft,A«))cos->  ,  (3.6) 

where  <f>  and  8  are  functions  of  the  known  angles,  Pi,  P2,  and  A 8  and  where  the  sgn  function 


is  defined  as  follows: 


f  —  1  if  x  <  0, 

sgn(x)  =  <  0  if  1  =  0,  and  (3.7) 

(  1  otherwise. 

Equation  (3.4)  may  also  be  used  to  find  sin#; 

.  ...  .  .  cos  P2  -  cos  Pi  cos  A9 

=  ,t„a»c0  (38) 

Equations  (3.5)  and  (3.6y  cannot  be  applied  if  sin  A#  =  0.  This  is  related  to  the  cases 
where  the  rotation  is  0°  and  180°.  Since  neither  of  these  cases  provides  an  effective  geometry 
change,  neither  of  these  cases  will  be  allowed. 

Furthermore,  if  cos<£  =  0,  equations  (3.6)  and  (3.8)  cannot  be  applied.  This  is  the  case 
of  the  source  being  directly  above  or  below  the  array.  If  <f>  =  90°,  then  both  Pi  =  90°  and 
P2  =  90°,  and  6  is  indeterminate.  Any  value  of  9  will  satisfy  equations  (3.3)  and  (3.4)  in  this 
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4  Extrema  of  Error  Intervals 


In  practice  the  observed  signal  arrival  angle  and  the  rotation  angle  of  the  horizontal  line 
array  are  usually  known  within  certain  error  bounds.  In  the  case  of  small  errors(in  the  signal 
arrival  angles  and  rotation  angle) ,  a  closed-form  approximation  for  the  errors  in  the  azimuthal 
and  vertical  arrival  angles  has  been  obtained  in  Reference  [I].  In  the  case  of  moderate  or  large 
errors,  however,  no  such  approximation  is  possible.  Instead,  results  for  the  latter  case  can  be 
obtained  conveniently  through  analytic  expressions  for  the  error  bounds  on  the  azimuthal  and 
vertical  arrival  angles,  for  arbitrary  error  bounds  on  the  input  parameters(i.e.  the  signal  arrival 
angle  and  the  rotation  angle).  This  section  focusses  on  the  derivation  of  these  analytic  bounds. 


4.1  Method  of  Solution  and  Definitions 

All  of  the  local  minima  and  maxima  of  6  and  <f>  in  the  observation  error  region  must  be 
found  in  order  to  find  the  global  minimum  and  maximum.  The  local  minima  and  maxima 
can  be  determined  by  finding  the  values  of  the  observed  angles  where  the  derivatives  of  6  or  <f> 
vanish.  The  values  of  9  or  <j>  at  such  points  may  be  compared  with  other  such  values  and  values 
on  the  boundary  of  the  region  to  determine  the  global  minimum  or  maximum. 

To  find  the  extrema,  the  first  partial  derivatives  must  be  taken.  To  simplify  the  partial 
derivatives,  each  variable  will  be  considered  to  have  radians  as  units.  The  departure  in  units  is 
for  this  section  only. 

Since  and  /32  are  only  used  in  cosine  form  and  since  <f>  is  only  used  in  cos2  form,  the 
following  transformations  can  be  used: 

Bi  =  cos  ,  (4.1) 

B2  =  cos  ft,  (4.2) 

*(Bi,B2,  A0)  =  cos2  4(01,02,  A0),  (4.3) 

and  0(B1,B2iA0)  =  *(0i,ft, A*).  (4.4) 

Because  of  the  0’s  only  appearing  in  cosine  form  and  <f>  only  appearing  in  cos2  form,  the  domain 
of  the  8' s  can  be  restricted  to  [0,  jr]  and  the  range  of  <f>  can  be  restricted  to  [0,  |],  This  ensures 
that  the  transformations  are  invertable. 

In  addition  to  the  above,  define  the  limits  on  the  error  intervals  as  follows: 


B l m ,n  =  minimum  possible  Bi  value, 
Bi_„  =  maximum  possible  Bi  value, 


and 


B2  . 

*min 

B, 

A0min 

0' 


9  min 

9  max 

4f 


$ min 
$max 


minimum  possible  B2  value, 
maximum  possible  B2  value, 
minimum  possible  A 9  value, 
maximum  possible  A 9  value, 
computed  azimuthal  arrival  angle, 
lower  bound  on  computed  value  of  O', 
upper  bound  on  computed  value  of  9', 
computed  vertical  arrival  angle, 
lower  bound  on  computed  value  of  <f>' , 
upper  bound  on  computed  value  of  <f>'. 


The  angle  A0  may  be  restricted  to  the  open  interval  (O.ir),  which  excludes  the  non- valid 
values  of  A0  =  0  and  A 9  =  tt.  Recall  that  these  angles  imply  that  no  effective  geometry  change 
has  occurred.  If  A0  is  in  (-jr,0),  it  may  be  modelled  by  using  |A0|  in  the  calculations  and 
negating  the  results  for  9'.  The  modification  of  the  values  may  be  represented  by  the  following 
algorithm. 

9'  :=  -9' 
t  :=  9min 
9  min  9ma  x 

9  min  t 

The  symbol  :=  represents  an  assignment  and  t  is  used  as  a  temporary  variable. 

Consider  3>(or  0)  as  a  mapping  from  a  three  dimensional  real  space  to  a  one  dimensional 
real  space.  The  domain  of  possible  values  for  Bj,  B2,  and  A 9  forms  a  cube  in  three  dimensions. 
The  global  minimum  and  maximum  of  $(or  0)  may  be  found  by  comparing  all  local  minima  and 
maxima  on  the  domain.  Global  minimum  and  maximum  refer  to  the  minimum  and  maximum 
over  the  entire  domain  while  the  local  minimum  and  maximum  refer  to  minimum  and  maximum 
on  an  arbitrarily  small  connected  subset  of  the  domain. 

All  such  local  minima  and  maxima  must  fall  either  in  the  interior  or  on  the  boundary  of 
the  cube.  A  minimum  or  maximum  may  occur  at  any  point  where  a  derivative  fails  to  exist. 
Otherwise,  if  a  point  in  the  interior  is  a  local  minimum  or  maximum,  all  directional  derivatives 
passing  through  that  point  must  vanish  there.  If  a  point  on  the  boundary  is  a  local  minimum 
or  maximum,  the  conditions  are  slightly  different. 

Let  x  be  a  point  on  the  boundary  where  a  local  minimum  or  maximum  occurs,  x  may 
be  a  vertex,  lie  on  an  edge  (but  not  a  vertex),  or  lie  on  a  face  of  the  cube  (excluding  edges 
and  vertices).  Suppose  x  lies  on  an  edge.  Two  observed  parameters  are  held  constant  and 
the  partial  derivative  with  respect  to  the  third  must  vanish  at  x.  Suppose  x  lies  on  a  face, 
but  not  on  an  edge.  In  this  case,  one  observed  parameter  is  held  constant  and  the  partial 
derivatives  with  respect  to  the  other  two  must  vanish  at  x.  These  conditions  are  necessary  but 
not  sufficient  for  the  existence  of  a  local  minimum  or  maximum. 

The  following  subsections  contain  the  analysis  showing  necessary  conditions  for  local  min¬ 
ima  and  maxima  of  0  and  $  to  exist.  This  is  done  by  finding  the  partial  derivatives  of  and 
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©  and  equating  them  to  zero.  Equating  each  of  the  derivatives  to  zero  implies  that  a  set  of 
conditions  must  be  satisfied.  Local  minima  and  maxima  may  be  found  by  using  various  com¬ 
binations  of  the  implied  conditions.  Finally,  all  of  the  local  minima  and  maxima  are  compared 
to  find  the  global  minimum  and  maximum. 

As  shown  in  the  following  sections,  the  interior  of  the  cube  may  contain  a  maximum  of  $ 
but  no  minimum.  The  maximum  may  be  in  the  interior  only  if  Pi  =  fa  =  f  can  be  satisfied 
there.  There  is  no  local  minimum  or  maximum  of  0  in  the  interior.  However,  if  the  cube 
contains  f}\  =  P2  =  §  1  then  0  is  indeterminate  there  and  hence  Sm*n  =  0  and  6maz  =  ?r. 

4.2  Local  Minima  and  Maxima  of  The  Vertical  Arrival  Angle 

As  shown  in  the  previous  section,  the  solution  for  <f>  may  be  written  as 


4>{PuPi,  AS)  =  cos-1 


COS  P\  COS  A  S  -  COS  Pi 

sin  AS 


^  +  COS2  Pi . 


Therefore,  $  may  be  expressed  as 


*(Bi,B2iAS)=  ( 


Bi  cos  AS  —  B2\2 


sin  AS 


)  +Bf. 


In  subsections  4.2.1  and  4.2.2,  the  location  of  the  minimum  and  maximum  of  <f>  (which  is 
the  same  as  of  <£)  will  be  determined. 

4.2.1  Local  Minima  and  Maxima  in  the  Interior  of  the  Domain 

Let  x  be  a  point  in  the  interior  of  the  domain  such  that  $  has  a  local  minimum  or  maximum 
at  x.  Then  each  directional  derivative  of  $  at  x  must  vanish  or  fail  to  exist.  In  particular,  the 
partial  derivatives  must  vanish  or  fail  to  exist. 

Each  of  the  partial  derivatives  can  be  found  as  follows.  The  partial  derivative  with  respect 
to  Bi  is 

nBi  -  B2  cos  AS 

8BT  “  2  »i„=  A9  '  (0) 


If  d$/dBi  =  0  then 


Bi  =  B2  cos  AS. 


The  partial  derivative  of  $  with  respect  to  B2  may  be  written  as 

3$  _  B2-B1COSAS 
3B2  sin2  AS 

If  d$/d B2  =  0  then  Bj,  B2,  and  AS  must  satisfy 

B2  =  Bi  cos  AS. 


(4.10) 


vv: 


Finally,  consider  the  derivative  of  $  with  respect  to  AO: 

3$  _ /BicosA^  -  B2\  /BjeosAS  -  Bi 


3A  0 


-0 


sin  Ad 


0(! 


sin2  Ad 


)• 


(4.11) 


If  3$/3Ad  =  0,  then  either 


B2  =  Bi  cos  Ad 


(4.12) 


or 


Bi  =  B2  cos  Ad. 


(4.13) 


Note  that  each  of  the  derivatives  exists  everywhere  on  the  domain. 

In  order  that  a  local  minimum  or  maximum  occur  at  x,  each  of  the  partial  derivatives 
must  vanish.  This  implies  that  equations  (4.8)  and  (4.10)  as  well  as  equations  (4.12)  or  (4.13) 
must  be  satisfied.  The  latter  two  equations  are  exactly  the  same  as  the  first  two  equations. 
In  order  that  (4.8)  and  (4.10)  be  satisfied,  both  Bi  and  B2  must  vanish  which  implies  that 
Pi  =  p2  =  j.  Hence  the  only  possible  local  maximum  or  minimum  in  the  interior  can  occur  at 
the  point  where  Pi  =  P 2  =  §•  This  corresponds  to  <f>  =  j  which  can  only  be  a  maximum  and 
therefore  a  local  minimum  may  only  occur  on  the  boundary. 


4.2.2  Local  Minima  and  Maxima  on  the  Boundary  of  the  Domain 


Possible  extrema  on  the  faces  of  the  domain  cube  will  now  be  considered.  The  faces  of  the 
cube  may  be  grouped  into  three  disjoint  sets  with  each  set  corresponding  to  faces  with  one  of 
the  three  parameters  being  held  constant.  The  parameter  held  constant  may  be  at  its  minimum 
or  maximum.  Each  set  is  discussed  separately  with  the  discussions  enumerated  1  through  3. 


1.  Consider  the  two  faces  where  Bi  is  held  constant:  that  is,  Bj  =  Bimm  or  Bj  = 

The  partial  derivatives  3$/3B2  and  3$/3A 0  must  vanish  at  the  point  where  an  extremum 
occurs.  As  shown  in  the  previous  discussion,  each  of  these  derivatives  vanishing  implies  that 
one  of  equations  (4.10),  (4.12),  or  (4.13)  holds.  d$/dA0  =  0  implies  that  either  equation  (4.12) 
or  equation  (4.13)  must  hold.  3$/3B2  =  0  implies  that  equation  (4.10)  must  be  satisfied  which 
is  the  same  as  equation  (4.12).  In  order  that  both  equation  (4.10)  and  equation  (4.13)  hold, 
Bj  =  B2  =  0  must  be  satisfied.  Hence  either  equation  (4.10)  or  Pi  —  P2  =  f  must  be  satisfied. 

2.  Repeat  the  same  procedure  with  B2  fixed:  that  is,  B2  =  B2m,„  or  B2  =  B2m„, . 
Therefore,  both  3$/3Bi  and  3<3>/3A 6  must  vanish  simultaneously  and  either  equation  (4.8)  or 
Bi  =  B2  =  0  [Pi  —  Pi  =  | )  must  be  satisfied. 

3.  Again,  repeat  the  same  procedure  with  AO  constant:  that  is,  3$/3Bi  =  3$/3B2  =  0. 
Equations  (4.8)  and  (4.10)  must  be  satisfied:  that  is,  Bi  =  B2  =  0  [Pi  =  P2  =  §)• 

The  conditions  for  local  extrema  on  the  edges  will  be  approached  in  a  similar  fashion. 
Each  edge  corresponds  to  two  observed  parameters  being  held  constant.  There  are  three  basic 
sets  of  edges  with  each  set  containing  four  edges;  each  set  refers  to  the  four  edges  corresponding 
to  two  parameters  being  held  constant  at  their  extrema.  Each  set  is  discussed  separately  with 
the  discussions  enumerated  1  through  3. 
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1.  Consider  the  edges  with  Bi  and  B2  fixed:  that  is,  Bj  =  Bimin  or  Bi  =  Bim„  and 
B2  =  B2min  or  B2  =  B2m#I.  Then,  for  a  point  to  be  a  local  minimum  or  maximum,  the 
partial  derivative  3$/3A 9  must  vanish  and  therefore  either  conditions  (4.12)  or  (4.13)  must 
be  satisfied. 


2.  Now  consider  Bi  and  A9  fixed  with  3$/3B2  vanishing.  3$/3B2  vanishing  implies  that 
equation  (4.10)  must  be  satisfied. 

3.  Finally,  consider  the  edges  with  B2  and  A 9  fixed.  3$/3Bi  vanishing  implies  that 
equation  (4.8)  must  be  satisfied. 

In  each  of  the  above  conditions  where  an  equation  (  (4.8),  (4.10),  (4.12),  or  (4.13) )  must  be 
satisfied,  only  one  point  satisfying  these  must  be  checked  since  $  is  constant  on  each  of  curves 
represented  by  these  equations.  Only  the  vertices  are  left  to  be  checked  for  possible  minima 
and  maxima.  Hence  only  a  finite  number  of  local  minima  and  maxima  must  be  compared  to 
find  the  global  minimum  and  maximum. 


4.3  Local  Minima  and  Maxima  of  Azimuthal  Arrival  Angle 


The  limits  of  the  error  interval  for  9  may  be  found  in  much  the  same  manner  as  with  <f>. 
The  azimuthal  angle  6  is  given  by 


0{Pi,P2,A0)  =  cos  1  ^ 


COS  fi\  ' 

cos <f>(fix, P2,  AQ), 


(4.14) 


Therefore,  0  is  given  by 


0(Bi,B2,  A0)  =  cos' 


(4.15) 


As  with  the  solution  for  the  bounds  of  the  error  interval  for  $,  the  partial  derivatives  of 
0  are  required.  The  partial  derivatives  of  0  are  given  by 


30  B2 

3BX  ~  $(Bi,B2,A0)sinA0’ 

30  _  _ Bi _ 

3B2  $(Bi,B2,  A0)sin  A6  ’ 

30  Bi(Bi  -  B2cos  A9) 

and  t  —  =  — — — — — — —5 


(4.16) 


(4.17) 


dA 9  $(Bi ,  B2 ,  A9)  sin2  A9 


(4.18) 


Equating  each  to  zero  implies  that  one  of  a  number  of  conditions  must  hold.  If  30/3Bj  =  0 


B2  =  0 


(4.19) 


which  implies  that  fit  =  f . 
If  30/3B2  =  0  then 


Bi  =  0 


(4.20) 


i 


© 


which  implies  that  8\  —  f  • 
If  80/ 8  AO  =  0  then 


Bi  (B2  cos  A0  -  Bi)  =  0. 


(4.21) 


Each  of  the  derivatives  fail  to  exist  if  $  =  0  or  sin  AO  =  0.  $  =  0  is  possible  only  if 
81  —  Pi  —  f  •  sin  AO  =  0  is  not  allowed. 

Again,  consider  solutions  for  the  extrema  on  the  vertices,  edges,  faces  and  interior  of  the 
domain  cube. 

4.3.1  Local  Minima  and  Maxima  in  the  Interior  of  the  Domain 

In  order  that  an  extreme  point  occur  in  the  interior  of  the  domain,  conditions  (4.19), 
(4.20),  and  (4.21)  must  be  satisfied  at  the  point.  This  would  imply  that  B2  =  B2  =  0  and 
therefore,  81  =  Pi  =  j  must  be  satisfied.  This  would  mean  that  $  =  0  and  therefore  0  would 
be  indeterminate. 

Hence,  0m,„  =  0  and  9max  =  *  if  81  =  Pi  —  §  i®  m  the  interior. 

4.3.2  Local  Minima  and  Maxima  on  the  Boundary  of  the  Domain 

As  in  subsection  4.2.2,  each  of  the  faces  and  edges  may  be  considered  separately  for  possible 
extrema.  There  are  three  sets  of  faces  to  be  considered;  each  set  has  a  particular  parameter 
held  at  its  minimum  or  maximum.  The  sets  are  discussed  separately  with  the  discussions  being 
enumerated  1  through  3. 

1.  Consider  the  faces  with  Bi  fixed:  that  is,  Bi  =  BimiB  or  Bj  =  Bim„.  In  order  that  a 
local  minimum  or  maximum  occurs  on  a  face  corresponding  to  Bi  fixed,  80/8 Bi  =  80/8 AO  = 
0  and  equations  (4.20)  and  (4.21)  must  be  satisfied.  Satisfying  these  implies  that  the  local 
minimum  or  maximum  can  occur  only  if  Bi  =  0;  i.e.  8\  —  j. 

2.  Now  consider  the  faces  with  B2  fixed:  that  is,  B2  =  B2ibi.b  or  B2  =  B2W„.  Then 
8Q/8Bi  =  80/8 AO  =  0  at  an  extremum.  Together,  conditions  (4.19)  and  (4.21)  imply  that 
Bi  =  B2  =  0  and  therefore  81  —  Pi  ~  \  must  be  satisfied.  Again  this  is  the  case  where  0  is 
indeterminate. 

3.  Finally,  consider  the  faces  with  AO  fixed:  that  is,  AO  =  A0m,n  or  AO  =  AO  mag.  In 
order  that  there  be  a  local  minima  or  maxima  on  the  face,  equations  (4.19)  and  (4.20)  must 
be  satisfied.  This  implies  that  Bi  =  B2  =  0  or  81  =  P2  =  § .  This  is  the  case  where  0  is 
indeterminate. 

The  conditions  necessary  for  local  minima  or  maxima  in  the  interior  and  on  the  faces  have 
been  considered.  Now  consider  the  edges  in  three  sets.  The  discussions  of  the  three  sets  of  four 
edges  is  enumerated  1  through  3.  Each  set  refers  to  the  edges  defined  by  two  parameters  held 
at  either  of  their  extreme  values. 

1.  Consider  the  edges  with  Bi  and  B2  fixed:  that  is,  Bj  =  Bimiit  or  Bi  =  Bimat  and 
B2  =  B2m.n  or  B2  =  B2m„.  In  order  that  there  be  an  extremum  on  one  of  the  edges,  equation 
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(4.21)  must  be  satisfied.  This  implies  that  either  Bi  =  0  or  B2  cos  A9  =  Bi  must  be  satisfied. 
0  is  constant  on  the  curves  represented  by  each  of  these  conditions  which  means  that  only  one 
point  need  be  found  for  each  of  the  conditions. 

2.  Now  consider  Bi  and  A 9  fixed:  that  is  Bi  =  Bimin  or  Bi  =  Bimai  and  A 9  =  A<?mtn  or 
A 9  =  A9max.  Equation  (4.20)  must  be  satisfied  at  an  extremum  which  implies  that  Pi  =  |.. 

3.  Finally,  consider  the  edges  with  B2  and  A 9  fixed:  that  is,  B2  =  B2mir>  or  B2  =  B2m„ 
and  A 9  —  A 0mm  or  A 9  =  A9maz.  The  existence  of  an  extreme  point  on  the  edges  implies  that 
equation  (4.19)  and  therefore  /?2  =  f  must  be  satisfied. 

The  globed  extremum  of  ©(and  9)  may  occur  at  any  point  in  the  domain  which  satisfies 
the  above  conditions  or  at  one  of  the  eight  vertices.  The  possible  points  must  be  compared  to 
find  the  global  minimum  and  maximum. 

4.4  Summary  on  Error  Bounds  for  <j>  and  9 

The  preceding  subsections  4.1,  4.2,  and  4.3  embrace  the  following  general  results: 

1.  The  domain  of  the  two  observed  signal  arrival  angles  (Bj  =  cos  Pi,  B2  =  cos  p2)  and 
the  rotation  angle(A0)  form  a  cube  in  three  dimensions.  The  global  minimum  and  maximum 
of  <j>  and  9  are  found  by  investigating  all  local  minima  and  maxima  over  the  entire  domain. 

2.  The  interior  of  the  cube  may  contain  a  global  maximum  but  no  local  minima  of  <f>. 

3.  It  is  possible  that  9  may  be  indeterminate  at  an  interior  point. 

4.  AH  other  local  minima  and  maxima  of  tj>  and  9  occur  on  the  boundaries  of  the  cube. 

5.  Finally,  only  the  local  extrema  at  a  finite  number  of  points  in  the  domain  need  be 
evaluated  and  compared  to  determine  the  global  minimum  and  maximum. 

The  detailed  descriptions  of  the  above  results  have  been  coded  in  Fortran  [2]  and  are  used 
in  the  next  section  to  obtain  graphs  of  maximum  error  bounds  for  the  azimuthal  angle  as  well 
as  upper  and  lower  error  bounds  for  the  vertical  arrival  angles. 
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5  Graphic  Example 


The  previous  section  presented  analytic  error  bounds  on  the  azimuthal  and  vertical  arrival 
angles  based  on  given  observation  error  bounds.  The  effect  of  observation  errors  on  the  com¬ 
puted  vertical  arrival  angle  <f>  is  considered  by  graphing  the  error  bounds  on  <f>  as  a  function  of 
the  azim  hal  angle  8.  Similarly,  the  effect  on  the  computed  azimuthal  angle  8  is  considered  by 
graphing  the  maximum  error  on  6  as  a  function  of  8.  Each  is  graphed  for  a  specified  rotation 
angle,  rotation  angle  error,  and  signal  arrival  angle  error.  Other  examples  are  presented  in 
Appendix  A. 

In  subsection  5.1,  the  errors  in  the  observed  data(the  signal  arrival  angles  and  the  rotation 
angle)  are  first  defined.  Next,  in  subsection  5.2,  the  definitions  of  the  the  bounds  on  the  input 
parameters^,  Pi,  and  Ad)  and  the  output  parameters  (<f>  and  8)  are  given.  Then,  in  subsection 
5.3,  the  results  for  the  case  of  small  observation  errors  are  presented.  Finally,  in  subsection  5.4, 
the  results  are  given  for  the  case  where  the  restriction  of  small  observation  errors  is  violated. 
For  consistency  with  the  results  in  Reference  [l],  all  angles  are  assumed  to  be  in  degrees. 

5.1  Errors  in  Observed  Data 

The  observation  errors  in  the  rotation  and  arrival  angles  may  be  defined  as  in  the  two 
following  subsections.  The  error  in  rotation  angles  is  constant  while  the  error  in  the  arrival 
angles  is  to  be  considered  constant  in  the  cosine  space. 

5.1.1  Errors  in  Observed  Signal  Arrival  Angles 

The  absolute  error  in  cos/9  is  assumed  to  be  independent  of  /?  and  is  denoted  by: 

e  =  sinc7o-  (5.1) 

With  the  definition  in  (5.1),  the  error  in  cos£  is  constant,  whereas  the  error  in  P  varies 
as  a  function  of  p.  Using  this  definition,  represents  the  error  in  P  when  p  is  equal  to 
90° (broadside).  Reference  [1]  shows  that  the  maximum  error  in  P  is  approximately  given  by 
I0.7y/oo.  This  occurs  near  the  endfire  direction. 

Therefore,  the  error  interval  for  cos  P  is  [Bmtn,  Bmaij  where 


and  B„ 


maz(-l,cos/?  -  e), 
m»n(l,cos/?  +  e). 


5.1.2  Errors  m  Rotation  Angle 


Further,  the  absolute  error  in  A8  is  denoted  by  h\  therefore,  the  error  interval  for  A 9  is 
given  by  [AO  -  h,  AO  -f  A]. 

5.2  Definitions  of  Bounds  on  Various  Parameters 


The  following  definitions  will  be  used: 


Pimin  =  cos-1  (mt'n((cos^i +sinff0),l)), 

Pima*  =  cos-1  (mai((cos^i  -  sinff0),-l)), 

02min  =  cos-1(m«'n((cos^2 +  sinff0),l)), 
p2max  =  cos-1  (mai((cos^2  -  sin«r0),-l)), 

A0mtn  =  A 0  -  h, 

Mmaz  =  A  e  +  h, 

O'  =  computed  azimuthal  arrival  angle, 

Omin  =  lower  bound  on  computed  value  of  8' , 

8 max  —  upper  bound  on  computed  value  of  O', 

4>  =  computed  vertical  arrival  angle, 

4>min  —  lower  bound  on  computed  value  of  <f>', 
and  <t>max  —  upper  bound  on  computed  value  of  <f>' 

where  ct0  is  assumed  to  lie  in  the  interval  [0,90).  Care  has  been  taken  through  the  use  of  the 
minimum  and  maximum  functions  to  ensure  that 


1  cos  01,,  J 

<  1, 

(5.4) 

l«>  spimj 

<  1, 

(5.5) 

1  COS  ^2m,nl 

<  1, 

(5.6) 

and  |cos£2mJ 

<  1. 

(5.7) 

The  maximum  error  in  0  is  defined  by 

terror  —  0 

,e'-0min). 

(5.8) 

5.3  Results  for  Small  Observation  Errors 

Reference  [l]  utilized  a  small  error  approximation  to  arrive  at  a  closed-form  solution  for 
the  errors  in  4>  and  0.  The  small  error  requirement  used  in  the  earlier  report  can  be  formally 
written  as: 

(5.9) 
15.101 


sin<7o  <?C  |sinA0|cos$ 

.  i  180,  .  A/11 

and  h  -  sinA0  . 


✓  / 


In  order  to  compare  the  results  obtained  by  the  method  using  the  extrema  of  error  intervals 
with  those  obtained  in  Reference  [l] (the  small  observation  error  approximation),  graphs  are 
presented  below  for  parameters  that  satisfy  Equations  (5.9)  and  (5.10).  This  comparison  serves 
to  validate  the  case  of  small  observation  errors  in  Reference  [l]. 

Figure  5.1  is  a  graph  of  the  maximum  error  in  the  estimated  azimuthal  angle  versus  the 
azimuthal  angle  for  A 6  =  20°,  ao  =  1°,  and  h  =  0°(From  the  comments  on  page  13,  a  value 
of  <ro  =  l°corresponds  to  an  error  in  fi  of  about  l°at  broadside  increasing  to  about  10.7°at 
the  endfire  directions).  As  the  vertical  angle  <f>  becomes  steeper,  the  effect  on  the  error  in 
the  computed  azimuthal  angle  becomes  larger  and  therefore  the  observation  errors  at  these 
higher  angles  have  a  correspondingly  larger  value.  Further,  the  largest  errors  occur  in  the 
endfire  directions,  where  cos  6  =  1,  whereas  the  minimum  errors  occur  near  broadside,  where 
cos/?  =  0.  Figures  A.l  to  A.6  in  Appendix  A  contain  similar  graphs  for  A8  =  45°,  90°,  and 
135°.  These  show  a  decrease  in  maximum  error  for  A 6  =  90°. 
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Figure  5.1:  Maximum  absolute  error  in  computed  azimuthal  angles  for  A 8  =  20°,  ao  =  1°,  and 
h  =  0°. 


Figure  5.2  shows  the  effect  of  cos/?  observation  errors  on  the  calculated  vertical  arrival 
angle  for  the  parameters  A 9  =  20°,  <7o  =  1°,  and  h  =  0°.  The  error  bounds  4>min  -  4>'  and 
<l>maz  ~  <t>'  are  plotted.  The  maximum  errors  occur  near  broadside  whereas  the  minimum  errors 
occur  near  the  endfire  directions. 

Figure  5.3  shows  the  effects  of  observation  errors  in  rotation  angle  on  the  computed  az¬ 
imuthal  angle  with  the  parameters  A 9  =  20°,  <ro  =  0°,  and  h  =  1°.  For  A6  =  20°,  the  maximum 
error  in  the  azimuthal  angle  is  independent  of  the  vertical  arrival  angle.  Figures  A.7  to  A.ll  in 
Appendix  A  contain  more  graphs  with  A 9  =  45°,  90°,  and  135°  which  show  the  same  trends. 
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Figure  5.2:  Error  bounds  on  computed  vertical  angles  for  A0  =  20°,  cr0  =  1°,  and  h  —  0° 
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Figure  5.3:  Maximum  absolute  error  in  computed  azimuthal  angles  for  A6  —  20°,  ao  =  0°,  and 
h  =  1°. 
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The  errors  in  computed  vertical  arrival  angles  are  much  less  influenced  by  errors  in  A 6 
than  errors  in  cos  /?.  Figure  5.4  shows  the  influence  of  errors  in  computed  vertical  arrival  angles 
for  A 6  =  20°,  <To  =  0°,  and  h  =  1°.  This  differs  from  the  effect  on  computed  azimuthal  angles  in 
that  the  errors  in  computed  vertical  arrival  angles  are  dependent  on  the  vertical  arrival  angles. 
The  corresponding  graphs  with  A 6  =  45° ,  90°,  and  135°  in  Appendix  A  also  support  this  result. 

ROTATION  ANGLE  A8-  20.00* 


AZIMUTHAL  ANGLE  0  IN  DEGREES 

Figure  5.4:  Error  bounds  on  computed  vertical  angles  for  A 9  =  20°,  ct0  =  0°,  and  h  =  1°. 


As  shown  in  Appendix  A,  the  effect  of  errors  in  cos  /?  on  the  computed  azimuthal  arrival 
angles  may  be  minimized.  Of  the  rotation  angles  used  in  Figure  5.1  and  in  Appendix  A,  a  90° 
rotation  would  minimize  such  effects.  The  effect  of  errors  in  A 8  on  the  computed  azimuthal 
arrival  angles  are  also  minimized  at  A 8  =  90°. 

The  consequence  of  observation  errors  on  the  computed  vertical  arrival  angle  show  similar 
trends.  In  fact,  the  effect  of  errors  in  cos  ,5  on  the  computed  vertical  arrival  angle  is  minimized 
at  A 9  =  90°.  Also,  the  minimum  effect  of  observation  errors  in  AS  on  the  computed  vertical 
arrival  angle  occurs  at  AS  =  90°. 

Each  of  the  above  sets  of  example  parameters  was  chosen  to  agree  with  previous  work 
[1].  A  comparison  of  the  above  graphs  as  well  as  those  in  Appendix  A  with  the  results  in 
Reference  [l]  showed  the  same  features  and  error  magnitudes.  The  error  magnitudes  are  the 
same  within  the  possible  accuracy  measured  from  the  graphs.  Hence,  the  method  of  extrema 
of  error  intervals  validates  the  results  using  the  small  observation  error  approximation. 
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5.4  Results  for  Moderate  Observation  Errors 

The  results  in  this  subsection  will  illustrate  the  point  that  the  approximate  error  bounds 
in  Reference  [1]  cannot  be  used  when  the  restrictions  in  Equations  (5.9)  and  (5.10)  are  violated. 
For  the  case  of  moderate  or  large  observation  errors,  the  method  using  the  extrema  of  the  error 
intervals  must  be  used. 

Four  graphs  showing  the  effect  of  moderate  observation  errors  have  been  prepared.  The 
effect  of  errors  in  the  observation  of  A 9  have  been  presented  for  A 9  =  20°.  Graphs  using  the 
results  obtained  in  previous  work[l]  have  been  included  for  comparison. 

The  error  parameters  h  =  4°  and  <tq  =  0°  violate  the  conditions  in  Equation  (5.10): 


180 

h  «  0.2 - j  sin  A9\. 

7T 


(5.11) 


The  values  of  these  parameters  are  used  to  illustrate  the  discrepancies  in  the  results  ob¬ 
tained  by  the  two  methods.  Figures  5.5  and  5.6  show  the  maximum  absolute  errors  in  the 
computed  azimuthal  arrival  angles;  Figure  5.5  makes  use  of  the  method  based  on  the  small  er¬ 
ror  approximation  whereas  Figure  5.6  uses  the  method  based  on  the  extrema  of  error  intervals. 
Similarly,  Figures  5.7  and  5.8  show  the  upper  and  lower  error  bounds  on  the  computed  vertical 
arrival  angles. 

The  graphs  on  the  following  two  pages  show  noticeable  differences  between  the  methods. 
The  approximate  method(for  small  observation  errors)  tends  to  underestimate  the  error  in 
azimuthal  angle  but  tends  to  overestimate  the  error  in  vertical  arrival  angle.  However,  it  can 
be  shown  that  for  other  values  of  the  parameters  h  and  Oq,  the  approximate  method  will  both 
underestimate  and  overestimate  the  error  in  azimuthal  angle  in  different  ranges  of  9.  Hence, 
the  above  results  indicate  that  the  method  of  extrema  of  error  intervals  developed  in  this 
paper  should  be  used  when  the  restriction  of  small  observation  errors  in  0  and  A 9  is  no  longer 
applicable. 
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6  Conclusion 


An  earlier  report  showed  that  if  a  Horizontal  Line  Array(HLA)  is  rotated  in  the  horizontal 
plane  to  a  new  orientation,  it  is  possible  to  solve  for  the  azimuthal  and  vertical  arrival  angles 
of  the  signal  [1],  The  same  report  gave  simple  expressions  for  r.m.s.  and  maximum  errors  in 
the  solution  when  the  observation  errors  are  small. 

In  this  report,  however,  the  theory  and  method  (of  extrema  of  error  intervals)  are  presented 
for  computing  error  bounds  on  the  solution  when  the  observation  errors  are  moderate  to  large. 
The  method  requires  evaluation  and  comparison  of  possible  solutions  at  a  finite  number  of 
points.  Though  maximum  errors  in  the  solution  are  usually  associated  with  maximum  errors  in 
the  input,  situations  arise  in  which  the  maximum  occurs  at  an  intermediate  value  of  the  input 
error.  This  method  has  been  implemented  in  Fortran  code [2], 

Graphs  of  the  errors  in  the  computed  azimuthal  and  vertical  arrival  angles  are  presented.  A 
comparison  has  been  made  with  the  previous  workfl]  for  the  cases  of  interest.  Good  agreement 
is  found  between  the  results  of  the  two  methods  for  small  observation  errors.  Discrepancies 
between  the  two  methods  are  shown  for  cases  where  the  restriction  of  small  observation  errors 
is  violated. 
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Appendix  A  Parametric  Analysis  of 
Maximum  Azimuthal  and  Vertical  Arrival 
Angle  Errors  for  Small  Observation  Errors. 


A  number  of  graphs  showing  the  effect  of  observation  errors  have  been  prepared  to  show 
the  effect  of  errors  in  cos  0  and  AO  observations.  These  are  presented  here. 

Graphs  have  been  presented  for  rotation  angles  AO  of  45°,  90°,  and  135°.  The  observation 
error  in  (3  given  by  cr0  is  defined  by 

cos (/3min)  =  mtn  ((cos/?  +  sin<70),  1) 
and  cos(/?maz)  =  max  ((cos/?  -  sin  <r0), -1) .  (A.l) 

The  observation  error  in  A 0,  given  by  h,  is  defined  by 

A0min  =  AO  -  h 

and  A0max  =  Ai 9  +  h.  (A. 2) 

For  each  rotation  angle,  graphs  are  presented  with  h  =  1°  or  ao  =  1°.  The  case  where 
both  and  h  are  non-zero  is  not  considered. 

Figures  A.l  through  A.6  show  the  maximum  absolute  errors  on  the  computed  azimuthal 
arrival  angles.  The  maximum  absolute  error  in  0  is  given  by 

0 error  =  rnax(Omat  -  O',  O'  -  0min)  (A. 3) 


where  O'  is  a  best  estimate  of  0.  Figures  A. 7  through  A. 12  show  the  upper  and  lower  error 
bounds  on  the  computed  vertical  arrival  angles  given  by  4>ma x  -  4>’  and  <j>'  -  where  <f>'  is  a 
best  estimate  of  the  vertical  arrival  angle  <f>. 
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Figure  A.l:  Maximum  absolute  error  in  computed  azimuthal  angles  for  A B  =  45°,  cto 
h  =  0°. 
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Figure  A.2:  Maximum  absolute  error  in  computed  azimuthal  angles  for  Ad  =  45°,  a0  —  0°,  and 
h  =  1°. 
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Figure  A.ll:  Error  bounds  on  computed  vertical  angles  for  A 6  =  135°,  a0  —  1°,  and  h  =  0°. 
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Figure  A. 12:  Error  bounds  on  computed  vertical  angles  for  A 9  =  135°,  =  0°,  and  h  =  1°. 
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